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1 .  INTRODUCTION 


When  a  planar  moving- incident  shock  wave 
impinges  on  a  sharp  compressive  corner  in  a  shock 
tube,  different  types  of  reflection  take  place 
depending  on  the  incident  shock  Mach  number  Ms  and 
the  corner  wedge  angle  0W.  The  four  possible  types 
of  reflection  in  pseudo-stationary  flows  are  illus¬ 
trated  in  Fig.  1:  (a)  regular  reflection  (RR) , 

(b)  single  Mach  reflection  (SMR) ,  and  in  Fig.  2: 

(a)  complex  Mach  reflection  (CMR),  and  (b)  double 
Mach  reflection  (DMR) .  A  RR  consists  of  a  two-shock 
configuration  (the  incident  shock  wave  I  and  the 
reflected  shock  wave  R) ,  as  shown  in  Fig.  1(a). 
However,  a  SMR  has  an  additional  third  shock  M, 
referred  to  as  a  Mach  stem,  and  a  slipstream  S, 
across  which  pressure  and  flow  direction  are  con¬ 
tinuous  but  all  other  flow  properties  are  discon¬ 
tinuous  [Fig.  1(b)].  The  confluence  point  T  of  the 
four  discontinuities  in  Fig.  1(b)  is  referred  to  as 
the  (first)  triple  point.  The  CMR  wave  system  is 
characterized  by  a  kinked  reflected  shock  wave 
[Fig.  2(a)],  whereas  the  DMR  has  two  reflected 
shock  waves  (R  and  R').  two  Mach  stems  (M  and  M') 
and  two  slipstreams  (S  and  S'),  as  shown  in  Fig. 
2(b). 

In  steady  supersonic  wind-tunnel  experiments 
only  RR  and  SMR  have  been  observed  and  this  distin¬ 
guishes  such  flows  from  nonstationary  one-dimen¬ 
sional  flows.  Ben-Dor  and  Glass  [Refs.  1,  2,  3] 
concluded  that  the  shock -wave  reflections  in  non¬ 
stationary  flows  also  depend  on  the  flow-deflection 
process  over  the  wedge  and  this  is  the  prime  reason 
for  the  existence  of  CMR  and  DMR.  The  incident 
shock  wave  is  reflected  by  the  wedge  surface, 
whereas  the  induced  nonstationary  flow  behind  it  is 
deflected  by  the  wedge  surface.  The  interaction 
between  these  two  processes,  which  take  place 
simultaneously,  causes  CMR  and  DMR  to  occur  in 
nonstationary  flows  in  addition  to  RR  and  SMR 
observed  in  steady  flows. 

Extensive  studies  on  the  transitions  between 
various  types  of  shock-wave  reflections  have  been 
done  by  many  researchers  in  the  past.  Griffith 
[Ref.  4]  has  given  a  recent  review.  Von  Neumann 
[Ref.  5]  studied  theoretically  the  RR  *■  MR  transi¬ 
tion  for  an  inviscid,  perfect-gas  flow  under  the 
concept  of  self-similarity.  He  deduced  the  well- 
known  detachment  criterion,  i.e.,  the  transition 
from  RR  to  MR  occurs  when  the  flow  deflection 
angle  0j  through  the  incident  shock  wave  exceeds 
the  maximum  flow  deflection  angle  02m  through  the 
reflected  shock  wave. 

Henderson  and  Lozzi  [Ref.  6]  proposed  the 
alternative  criterion,  which  is  referred  to  as  the 
mechanical-equilibrium  criterion  for  the  RR  ?  MR 
transition.  They  investigated  the  RR  t  MR  transi¬ 
tion  problem  by  using  both  a  wind  tunnel  and  a 
shock  tube.  This  criterion  is  attributed  to  the 
consideration  that  the  pressure  change  in  the 
system  should  be  gradual  at  the  transition  from 
RR  to  MR;  in  other  words,  the  system  remains  in 
mechanical  equilibrium  during  the  transition. 
Recently.  Hornung,  Oe’-  »l  and  Sandeman  [Ref.  7] 
tried  to  ♦’•Muce  .<w  criterion  with  a  length 
scale  in  tne  «  ste...  However,  their  criterion 
and  the  lengzr.  scale  have  r.ot  been  substantiated. 
According  to  the  sonic  criterion  proposed  by 
Horaung  et  al  [Ret.  7],  the  RR  £  MR  transition 
takes  place  at  the  sonic  deflection  angle,  the 


angle  of  deflection  at  which  the  flow  behind  the 
reflected  shock  wave  is  just  sonic  relative  to  the 
reflection  point.  Note  that  the  sonic  deflection 
angle  and  the  maximum  deflection  angle  are  very 
close  [see  Ref.  8].  Therefore,  the  transition 
lines  predicted  by  the  detachment  and  the  sonic 
criteria  are  almost  identical  and  ar^  difficult 
to  resolve  experimentally. 

The  experimental  data  of  Ben-Dor  [Ref.  1], 
Henderson  and  Lozzi  [Ref.  6],  Hornung  et  a)  [Ref. 

7],  Smith  [Ref.  9],  White  [Ref.  10]  and  others 
support  the  mechanical  equilibrium  criterion  as 
the  correct  RR  5  MR  transition  criterion  for 
steady  flows;  meanwhile  the  detachment  criterion 
is  accepted  as  the  correct  transition  criterion 
for  the  nonstationary  case.  Smith  [Ref.  9]  per¬ 
formed  intensive  shock-tube  experiments  for  shock- 
wave  reflections.  Bleakney  and  Taub  [Ref.  Ill 
revised  Smith's  data  and  mapped  out  the  region  of 
RR  (two-shock  pattern)  and  Mach  reflection  MR 
(three-shock  pattern) .  Some  disagreement  between 
the  detachment  criterion  and  Smith's  experimental 
results  were  found,  i.e.,  RR  persisted  beyond  the 
limit  determined  by  the  detachment  criterion  at 
low  incident  shock -wave  Mach  numbers.  This  per¬ 
sistence  of  RR  beyond  the  theoretical  limit  was 
observed  by  Henderson  and  Siegenthaler  [Ref.  12] 
at  similar  conditions  and  a  cross-over  at  a  higher 
pressure  ratio. 

The  attempt  of  Henderson  and  Lozzi  [Ref.  6]  to 
substantiate  the  mechanical  equilibrium  criterion 
revealed  a  remarkable  anomaly  between  their  results 
from  wind-tunnel  and  shock -tube  experiments.  They 
found  that  RR  continued  to  exist  in  the  shock-tube 
experiments  beyond  the  limit  predicted  by  both  the 
detachment  and  the  mechanical  equilibrium  criteria. 
They  resolved  the  anomaly  by  postulating  that  the 
RR  configurations  observed  beyond  the  theoretical 
limit  were  undeveloped  DMR  configurations  in  which 
all  the  shock  waves,  slipstreams  and  triple-points 
typical  of  a  well-developed  DMR  were  too  close 
together  to  be  observed.  Ben-Dor  and  Glass  [Rers. 

2,  3]  and  Ando  [Ref.  13]  criticized  this  hypothesis, 
based  on  the  experimental  results,  as  doubtful. 

If  a  moving- incident  shock  wave  encounters 
concave  or  convex  wedges  instead  of  a  sharp  wedge, 
a  transient  effect  arises  in  the  RR  X  MR  transition 
phenomenon.  Since  the  wedge  angle  0W  varies  along 
concave  and  convex  wedges,  then  the  entire  process 
will  be  transient  and  the  flow  is  truly  nonstation¬ 
ary.  The  existing  criteria  (detachment  and  mechan¬ 
ical  equilibrium)  for  the  RR  *■  MR  transition  cannot 
be  applied  to  the  truly  nonstationary  case,  since 
they  were  based  on  stationary  flow.  Ben-Dor  and 
Glass  [Ref.  2]  hypothesized  that  in  truly  nonsta¬ 
tionary  flows  there  existed  a  hysteresis  between 
the  RR  ■*  MR  and  the  MR  -*•  RR  transitions.  Ben-Dor 
et  al  [Ref.  14]  and  Itoh  et  al  [Ref.  15]  investi¬ 
gated  the  RR  ^  MR  transition  in  truly  nonstationary 
flows  and  confirmed  this  hypothesis. 

When  RR  terminates,  three  different  types  of 
Mach  reflection  can  occur  in  pseudo-stationary 
flows.  They  are  CMR  and  DMR,  in  addition  to  SMR. 

As  noted  previously,  CMR  first  noticed  by  Smith 
[Ref.  9],  and  DMR  first  observed  by  White  [Ref.  10], 
take  place  as  a  result  of  the  interaction  between 
the  reflection  process  at  the  wedge  surface  and 
the  flow-deflection  process  over  the  wedge  corner, 
which  is  peculiar  to  time-dependent  flows  in  a 


shock  tube.  White  [Ref.  10]  investigated  the  addi¬ 
tional  reflection  patterns  (CMR  and  DMR)  and  noticed 
that  both  of  them  were  associated  with  the  flow 
behind  the  reflected  shock  wave  which  was  supersonic 
relative  to  the  triple-point  T  (i.e.,  M2T  >  1.0). 

A  CMR  is  characterized  by  a  straight  reflected 
shock  wave  at  the  triple-point  T  that  is  connected 
to  a  curved  bow  wave  through  a  smooth  or  sharper 
kink  K  [Fig.  2(a)],  This  pattern  is  distinguished 
from  that  of  SMR  which  has  no  kink  in  the  reflected 
shock  wave  over  its  entire  region  [Fig.  1(b)].  A 
DMR  also  has  a  straight  reflected  shock  wave  at  T, 
and  a  secondary  reflected  shock  wave  R'  emanating 
from  a  sharp  kink  T',  a  secondary  Mach  stem  M' , 
and  a  secondary  slipstream  S'  originating  at  the 
second-triple-point  T'  [Fig.  2(b)].  Henderson  and 
Lozzi  [Ref.  6]  suggested  that  a  band  of  compression 
waves  must  exist  in  a  CMR.  These  compression  waves 
finally  converge  to  a  shock  wave  and  DMR  occurs 
when  the  flow  behind  the  reflected  shock  wave  with 
respect  to  the  first-triple-point  T  is  supersonic. 
Gvozdeva  et  al  [Ref.  16,  17]  investigated  the 
mechanism  for  the  formation  of  DMR  and  erroneously 
suggested  that  DMR  occurs  as  a  result  of  the 
excitation  of  internal  degrees  of  freedom  of  the 
gas  molecules  behind  the  reflected  shock  wave  and 
the  Mach  stem.  Bazhenova  et  al  [Ref.  18]  obtained 
the  transition  boundary  lines  for  SMR,  CMR  and  DMR 
experimentally.  However,  their  experimental  data 
were  limited  to  lower  Mach  number  and  did  not  cover 
a  wide  range  of  interest . 

Ben-Dor  [Ref.  1]  and  Ben-Dor  and  Glass  [Refs. 

2,  3]  established  the  transition  criteria  which  were 
initiated  by  Law  and  Glass  [Ref.  19]  for  the 
SMR  £  CMR  and  the  CMR  t  DMR  transitions.  They 
suggested  that  the  SMR  2  CMR  transition  occurred 
when  the  flow  behind  the  reflected  shock  wave  R  was 
sonic  with  respect  to  the  first-triple-point  T 
(i.e.,  M2t  =  1.0)  and  that  the  CMR  t  DMR  transition 
took  place  when  the  flow  behind  the  reflected 
shock  wave  R  was  sonic  relative  to  the  kink  K  (i.e., 
M2k  =  1.0).  It  must  be  noted  that  the  kink  K,  or 
the  second-triple-point  T',  is  formed  as  a  result 
of  the  interaction  between  the  shock  reflection  and 
the  flow-deflection  processes.  Since  the  above 
mentioned  criteria  for  the  SMR  t  CMR  and  the 
CMR  t  DMR  transitions  are  deduced  without  consider¬ 
ing  that  interaction  process  explicitly,  their 
physical  interpretations  are  still  somewhat  ambigu¬ 
ous,  even  though  the  transition  lines  predicted 
by  these  criteria  agree  quite  well  with  experimental 
data  for  some  gases. 

Ben-Dor  [Ref.  1]  and  Ben-Dor  and  Glass  [Refs. 

2,  3]  mapped  out  analytically  the  regions  of  RR, 

SMR,  CMR  and  DMR  in  the  (Ms-94) -plane,  where  Ms  is 
the  incident  shock  Mach  number  and  9^  the  effective 
wedge  angle  (9i  =  9W  +  x>  where  x  is  the  triple¬ 
point-trajectory  angle).  They  conducted  the 
experiments  for  the  various  types  of  reflection  in 
nitrogen  (as  a  diatomic  gas)  and  argon  (as  a  mon¬ 
atomic  gas)  over  the  range  of  Ms  from  1.9  to  8.0. 
Following  the  work  of  Ben-Dor  and  Glass,  Ando  [Ref. 
13]  delineated  the  regions  of  the  various  types  of 
reflection  for  carbon  dioxide  (as  a  linear  tri- 
atomic  gas)  in  the  (Ms-9w) -plane  and  performed  the 
experiments  over  the  range  of  Ms  from  1.8  to  10.2. 

It  is  natural  to  consider  that  the  real-gas 
effects  due  to  the  excitation  of  the  internal 
degrees  of  freedom,  vibration,  dissociation  and 
ionization  may  influence  the  shock-wave  system  in 


case  of  strong  shock  waves.  Gvozdeva  et  al  [Ref. 

16]  studied  Mach  reflection  in  various  perfect 
gases  with  different  specific-heat  ratios  y,  and  in 
vibrational ly  excited  carbon-dioxide  and  air.  They 
noted  that  an  increase  of  the  incident  shock  velo¬ 
city  led  to  a  decrease  of  the  reflection  angle  u>' 
[Figs.  1,  2],  and  becoming  negative  under  certain 
conditions.  Hornung  et  al  [Ref.  7]  investigated 
the  RR  *  MR  transition  for  dissociating  nitrogen 
and  carbon  dioxide,  and  ionizing  argon.  In  addi¬ 
tion  to  real-gas  effects,  they  gave  an  explanation 
for  the  persistence  of  RR  beyond  the  limit  predicted 
“by  the  detachment  criterion  in  terms  of  the  displace 
ment  thickness  of  the  boundary  layer  on  the  wall 
surface  in  pseudo-stationary  flow. 

Ben-Dor  [Ref.  1]  and  Ben-Dor  and  Glass  [Refs. 

2,  3]  studied  analytically  and  experimentally  the 
domains  and  boundaries  between  various  types  of 
reflection  patterns  in  dissociating  nitrogen  and 
ionizing  argcn  under  equilibrium  conditions.  They 
concluded  that  real -gas  effects  had  a  significant 
influence  on  the  size  of  the  regions  and  their 
boundaries;  hence  the  shock -wave  pattern  depended 
on  the  initial  thermodynamic  state  of  the  gas 
(i.e.,  temperature  To  and  pressure  P0)  in  addition 
to  Ms  and  9w.  They  conducted  shock-tube  experi¬ 
ments  covering  a  wide  range  of  Ms  and  9W;  however 
these  experiments  were  not  directed  sufficiently 
towards  the  transition  lines  of  the  various  types 
of  reflection.  Ando  [Ref.  13]  calculated  the 
regions  corresponding  to  the  four  types  of  reflec¬ 
tion  in  carbon  dioxide  using  five  models  of  the 
equation  of  state  from  complete  dissociational- 
vibrational  equilibrium  to  a  perfect  (frozen)  gas. 
Based  on  his  experimental  results  in  the  (Ms-9w)- 
plane,  he  concluded  that  real-gas  effects  played 
no  part  in  predicting  the  boundaries  between 
various  types  of  reflection,  and  hence  all  transi¬ 
tion  boundaries  had  to  be  determined  on  the  basis 
of  a  perfect  gas.  Ando  and  Glass  [Ref.  20]  hypo¬ 
thesized  to  explain  this  as  follows:  The  shock 
fronts  in  the  various  wave  systems  instantly 
responded  to  the  existing  state  of  the  gas.  Any 
further  vibrational  or  dissociational  excitation 
required  appropriate  relaxation  times  and  if  large 
would  not  be  a  contributing  factor  to  the  immediate 
shock  shapes  at  the  point  of  reflection  P  or  triple¬ 
point  T.  Their  explanation  appears  reasonable  for 
nitrogen  or  even  air  where  the  relaxation  lengths 
are  large.  However,  for  carbon  dioxide  this  is  not 
the  case.  Yet,  the  best  agreement  in  the  (Ms-9w)~ 
plane  occurred  for  C02  as  a  perfect  gas.  This 
casts  some  doubt  on  the  accuracy  of  the  present 
criteria,  even  though  they  are  the  best  available 
to  date.  This  will  be  discussed  further  by 
Shirouzu  and  Glass  [Ref.  36]  in  the  near  future. 

Note  that  if  the  relaxation  length  (i.e.,  the 
distance  behind  the  shock  wave  in  which  the  gas 
reaches  equilibrium  and  depends  on  the  gas  and  the 
initial  thermodynamic  state)  is  small  compared  with 
the  available  flow  length,  the  flow  is  in  equili¬ 
brium.  On  the  other  hand,  if  the  relaxation  length 
is  large  compared  with  the  region  of  interest,  the 
flow  is  frozen.  If  the  relaxation  length  is  com¬ 
parable  to  the  available  flow  length,  the  flow  will 
be  in  nonequilibrium.  Consequently,  real-gas 
effects  could  affect  the  location  of  K  and  T',  or 
R  and  R',  and  awaits  analytical  and  experimental 
verification. 

Several  numerical  techniques  have  been  proposed 
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to  solve  and  numerically  simulate  the  full  flow 
fields  involving  the  multiple  shock-wave  configura¬ 
tions.  Auld  and  Bird  [Ref.  21],  Schneyer  [Ref.  22], 
(Cutler  and  Shankar  [Ref.  23]  and  Shankar  et  al 
[Ref.  24]  performed  numerical  studies  on  regular 
and  single-Mach  reflections  using  different  numer¬ 
ical  techniques.  Recently  Booen  and  Needham  [Ref. 
25],  Book  et  al  [Ref.  26],  Glaz  [Ref.  27]  and 
Champney  et  al  [Ref.  28]  presented  results  of  num¬ 
erical  calculations  simulating  complex  and  double- 
Mach  reflections.  The  earlier  calculated  results 
[Ref.  22-24]  showed  satisfactory  agreement  with 
the  actual  experimental  results  in  the  overall  wave 
shapes  and  systems.  However,  they  gave  rather 
poor  predictions  for  the  flow  properties  such  as 
the  isopycnics  and  fair  agreement  for  the  density 
profiles  along  the  wall  [see  Ben-Dor  and  Glass, 

Ref.  29],  The  later  results  [Refs.  25-28]  show 
much  better  agreement  for  the  density  along  the 
wall,  and  improved  though  still  not  satisfactory 
isopycnic  shapes  and  locations.  Even  though  inten¬ 
sive  efforts  are  required  to  improve  the  numerical 
techniques  for  obtaining  solutions  with  acceptable 
costs,  the  numerical  methods  which  are  presently 
evolving  will  be  valuable  assets  in  predicting 
such  complex  shock-wave  flows  in  the  near  future. 

The  objectives  of  the  present  study  were  as 
follows : 

(1)  To  predict  the  domains  and  transition  boundaries 
of  the  various  types  of  reflection  for  perfect 
and  imperfect  air  in  thermodynamic  equilibrium 
in  the  (MS-0W) -plane. 

(2)  To  clarify  the  effects  of  the  specific-heat 
ratio  y  on  the  transition  boundaries  and  the 
shock-wave  configurations. 


of  the  confluence  point  [the  reflection  point  P  for 
RR,  Fig.  1(a),  or  the  triple-point  T  for  MR,  Fig. 
1(b)],  the  calculations  for  obtaining  the  flow 
properties  of  each  region  only  require  the  use  of 
oblique-shock-wave  relations  for  each  shock  wave. 
This  is  based  on  the  assumption  that  in  the  angular 
regions  around  the  confluence  point  the  flow  is  in 
a  uniform  state.  The  necessary  boundary  condition 
for  regular  reflection  (two-shock  theory)  is  that 
the  flow  deflections  across  the  incident  and  re¬ 
flected  shock  waves  must  be  equal  and  opposite, 
i.e. , 

Oj  +  S,  =  0  (1) 

The  boundary  conditions  for  Mach  reflection  (three- 
shock  theory)  are  that  the  total  flow  deflection 
across  the  incident  and  reflected  shock  waves  must 
be  equal  to  the  flow  deflection  across  the  Mach 
stem,  and  that  the  pressures  across  the  slipstream 
are  equal ,  i.e.. 


P2  =  P3  (3) 

In  a  MR,  it  is  necessary  to  give  a  value  to  the 
triple-point-trajectory  angle  x  as  an  initial  con¬ 
dition  in  order  that  the  system  of  equations  be 
closed  [Ben-Dor,  Ref.  1].  Law  and  Glass  [Ref.  19] 
proposed  an  empirical  method  for  predicting  the 
value  of  x  based  on  experimental  observations. 

They  assumed  a  straight  Mach  stem  normal  to  the 
wedge  surface  and  introduced  an  additional  inde¬ 
pendent  geometrical  relation  as: 


(3)  To  investigate  the  validity  and  the  restrictions 
of  existing  transition  criteria  for  the  various 
types  of  reflection. 

(4)  To  estimate  the  contribution  of  real-gas  effects 
to  the  shock-wave  configurations  by  comparing 
the  analytical  results  with  the  available  ex¬ 
perimental  data. 


ANALYTICAL  CONSIDERATION 


2.1  Analysis  for  Pseudo-Stationary  Flow 


The  above  assumption  is  good  in  the  range  5°  i  0W  S 
45°  [Ben-Dor,  Ref.  1].  However,  it  is  not  valid 
for  Ms  <  2  and  0W  £  5°.  This  restriction  will  be 
explained  subsequently. 

Since  the  analytical  formulations  of  the  two- 
shock  and  three-shock  theories  and  the  calculation 
procedures  are  well  known  [see  Ben-Dor,  Ref.  1,  and 
Law  and  Glass,  Ref.  19],  none  of  the  details  will 
be  given  here. 

Calculations  were  performed  for  the  following 
two  cases: 


Since  the  planar  incident  shock  wave  I  [Figs. 

1  and  2]  moves  with  constant  velocity  along  the 
wedge  surface,  the  entire  reflection  phenomenon  can 
be  considered  to  be  pseudo-stationary  in  a  frame 
of  reference  attached  to  the  confluence  point  of 
the  shock  waves.  That  is,  instead  of  three  inde¬ 
pendent  variables  x,  y  and  t,  the  phenomenon  may 
be  described  in  terms  of  x/t  and  y/t,  and  the  flow 
is  now  self-similar  [Bleakney  and  Taub,  Ref.  11], 
The  variables  x  and  y  may  be  measured  relative  to 
any  point  moving  with  constant  velocity  with  re¬ 
spect  to  the  wedge  corner,  and  time  t  is  measured 
from  the  instant  the  incident  shock  wave  has  passed 
through  the  point . 

In  a  pseudo-stationary  flow,  the  entire  shock- 
wave  configuration  remains  similar  and  grows 
linearly  with  time  from  the  instant  the  incident 
shock  wave  collides  with  the  wedge  corner.  If  we 
pay  attention  to  the  flow  fields  in  the  vicinity 


(1)  Perfect  air  with  a  specific-heat  ratio  y  =  7/5. 

(2)  Thermodynamic  equilibrium  imperfect  air,  taking 
into  account  rotation-vibration  coupling,  vibra¬ 
tional  excitation,  dissociation,  electronic  excita¬ 
tion  and  ionization.  In  this  case, the  equilibrium 
calculation  of  air  was  made  by  using  the  equilibrium - 
constant  technique  of  Hansen  [Ref.  30]  with  7 
species  (N2,  O2,  N,  0,  N+,  0+,  e~).  The  mole 
fractions  of  other  species  except  nitric-oxide  NO 
are  of  the  order  of  0.1  percent,  or  less  in  high- 
temperature  air  [see  Gilmore,  Ref.  31,  and  Laird 

and  Heron,  Ref.  32],  The  mole  fraction  of  NO  may 
become  as  much  as  10  percent  at  a  pressure  of  10 
atmospheres  and  a  temperature  of  about  5000  K. 
However,  this  nitric-oxide  does  not  strongly  in¬ 
fluence  the  resulting  thermodynamic  properties  of 
air.  Therefore,  NO  was  neglected  in  the  present 
equilibrium  calculation  for  simplicity.  The  mole 
fractions  of  molecular  oxygen  and  nitrogen  in  air 
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at  room  temperature  were  assumed  to  be  0.2  and  0.8, 
respectively.  For  the  imperfect-air  case,  the  gas 
is  considered  to  reach  a  new  equilibrium  state  imme¬ 
diately  behind  each  shock  wave  around  a  confluence 
point. 


Formation  and  Termination  Criteria  for  the 
Various  Types  of  Reflection  — 


In  order  to  predict  the  transition  boundaries 
and  the  regions  of  the  various  types  of  shock-wave 
reflection,  the  following  transition  criteria  were 
utilized  as  the  most  reasonable  criteria  (confirmed 
by  experiments)  to  date. 


The  detachment  criterion  proposed  by  von  Neumann 
[Ref.  5]  was  used  for  the  termination  criterion  of 
RR.  This  criterion  is  attributed  to  the  concept 
that  the  transition  from  RR  to  MR  occurs  when  the 
flow-deflection  angle  01,  through  the  incident 
shock  wave  exceeds  the  maximum  flow  deflection 
angle  02m,  through  the  reflected  shock  wave.  Ana¬ 
lytically,  this  can  be  expressed  as 


transition,  because  the  first  triple-point  T  and 
the  kink  K  (or  the  second  triple-point  T ’ )  coin¬ 
cided  at  the  RR  Z  MR  transition  boundary.  Conse¬ 
quently,  the  SMR  t  CMR  transition  line  and  the 
CMR  Z  DMR  transition  line  must  merge  at  the  RR  t  MR 
transition  boundary.  Further  analysis  for  the  pre¬ 
diction  of  the  location  of  K  (or  T')  at  higher 
wedge  angle  0W  is  still  required  in  order  to 
analytically  modify  the  CMR  £  DMR  transition  line  C 
[Fig.  3(a)], 


3 .  RESULTS  AND  DISCUSSIONS 

3.1  Analytical  Results  and  Discussions 


The  transition  boundary  lines  between  various 
types  of  nonstationary  oblique-shock -wave  reflec¬ 
tion  were  calculated  based  on  the  two-shock  and 
three-shock  theories.  The  aforementioned  formation 
and  termination  criteria  for  RR,  SMR,  CMR  and  DMR 
were  used  along  with  the  basic  assumption  for  the 
shock -wave  configurations. 


For  the  SMR  Z  CMR  and  the  CMR  Z  DMR  transitions, 
the  criteria  established  by  Ben-Dor  [Ref.  1]  and 
Ben-Dor  and  Glass  [Refs.  2,  3]  were  applied.  They 
suggested  that  the  SMR  £  CMR  transition  occurred 
when  the  flow  behind  the  reflected  shock  wave  R  was 
sonic  with  respect  to  the  first -triple-point  T,  and 
that  the  CMR  £  DMR  transition  took  place  when  the 
flow  behind  the  reflected  shock  wave  R  was  sonic 
with  respect  to  the  kink  K.  Therefore  the  SMR  £  CMR 
transition  criterion  is  expressed  as 

M2T  =  1.0  (6) 

and  the  CMR  Z  DMR  transition  criterion  as 


M2|(  =  1-0  (7) 

In  addition,  it  is  necessary  to  predict  the 
location  of  the  kink  K,  in  order  to  obtain  the  flow 
velocity  behind  the  reflected  shock  wave  with  res¬ 
pect  to  the  kink.  Law  and  Glass  [Ref.  19]  assumed 
that  the  kink  K  moved  with  the  same  horizontal 
velocity  as  the  induced  flow  behind  the  incident 
shock  I.  As  a  result,  the  ratio  L/L,  where  L  is 
the  horizontal  distance  between  the  first-triple¬ 
point  T  and  the  kink  K  (or  the  second  triple-point 
T'),  and  L  is  the  distance  between  T  and  the  vertex 
of  the  wedge  [see  Fig.  2],  becomes  equal  to  the 
density  ratio  across  the  incident  shock  wave  irre¬ 
spective  of  the  wedge  angle  0W,  that  is 


where, p0/pj  is  the  density  ratio  across  the  incident 
shock  wave  I. 


Bazhenova  et  al  [Ref.  18]  investigated  the 
ratio  L/L  experimentally  in  argon,  air  and  carbon 
dioxide  and  found  that  the  relation  (8)  was  valid 
only  in  the  range  8w  <  40°.  Their  experimental  data 
for  8W  >  40°  showed  that  the  ratio  L/L  decreased 
and  approached  zero  as  8w  increased.  This  fact  was 
confirmed  experimentally  by  Ando  and  Glass  [Ref.  20]. 
They  concluded  that  the  ratio  L/L  vanished  at  the 
critical  wedge  angle  corresponding  to  the  RR  Z  MR 


Figures  3(a)  and  3(b)  show  the  domains  and 
transition  boundaries  for  various  types  of  reflec¬ 
tion  in  the  (Ms -0^) -plane.  Note  that  the  effective 
wedge  angle  0^  is  the  sum  of  the  actual  wedge  angle 
8W  and  the  first  triple-point-trajectory  angle  x 
[see  Figs.  1,  2],  In  the  following,  the  solid 
boundary  lines  are  for  perfect  air  with  y  =  7/5, 
while  the  dashed  lines  are  for  air  in  thermodynamic 
equilibrium  (rotation-vibration  coupling,  vibrational 
excitation,  dissociation,  electronic  excitation  and 
ionization)  at  the  given  initial  pressure  P0  =  IS 
torr  and  temperature  T0  =  ZOO  K.  The  initial  values 
were  chosen  to  correspond  to  previous  work  and  ex¬ 
perimental  data.  Lines  A  and  B  show  the  RR  £  MR 
transition  lines  (8^  +  02m  -  0)  for  perfect  and  im¬ 
perfect  air,  respectively.  Lines  C  and  D  are  the 
CMR  *  DMR  transition  line  (M27  =  1.0),  while  lines 
E  and  F  are  the  SMR  £  CMR  transition  lines  (M2y  = 
1.0).  Lines  G  and  H  show  the  lines  of  Mi  =  1.0, 
below  which  no  reflection  (NR)  can  occur.  It  is 
clear  from  these  figures  that  the  boundary  lines  for 
perfect  air  (lines  A,  C,  E  and  G)  are  monotonic  and 
smooth.  As  the  incident  shock  Mach  number  Ms  in¬ 
creases,  these  lines  level  out  and  become  almost 
independent  of  Ms-  This  will  be  explained  subse¬ 
quently. 

It  is  evident  from  Figs.  3(a)  and  3(b)  that 
real-gas  effects  shift  every  boundary  line  down¬ 
wards  in  the  (Ms-8^) -plane.  That  is,  the  imperfect 
boundaries  begin  to  drop  from  their  perfect-gas 
values  even  at  very  low  incident  shock -wave  Mach 
number  Ms,  owing  to  vibrational  excitation.  The 
line  B,  for  instance,  starts  to  separate  from  line 
A  at  Ms  =  1.40  and  decreases  as  Ms  increases.  It 
must  be  noted  that  real -gas  effects  for  air  appear 
more  clearly  than  for  pure  nitrogen  because  of  the 
presence  of  about  20%  oxygen,  whose  internal  energy 
is  excited  at  lower  shock  Mach  numbers  than  that  of 
nitrogen.  Dissociation  also  occurs  at  lower  shock 
Mach  number  in  oxygen  compared  to  nitrogen.  At 
higher  values  of  Ms,  dissociation  of  oxygen  in  air 
occurs  first,  followed  by  nitrogen.  Electronic 
excitation  and  ionization  of  dissociated  atoms 
becomes  significant  at  still  higher  shock  Mach 
numbers.  Due  to  these  real-gas  effects,  the  im- 
peifect  lines  (B,  D,  F  and  H)  in  air  have  quite 
different  features  from  those  of  the  perfect-gas 
lines  (A,  C,  E  and  G) .  The  imperfect-gas  lines 
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not  only  decrease  but  also,  in  some  ranges,  increase 
with  increasing  Ms,  i.e.,  the  lines  undulate  as  Ms 
increases .  The  general  tendency  of  the  imperfect 
boundary  lines  can  be  understood  by  considering  the 
density  ratio  across  each  shock  wave.  The  density 
ratio  across  a  normal  shock  wave  in  imperfect  air 
increases  rapidly  but  undulates  within  certain 
ranges  of  shock  Mach  number.  This  trend  becomes 
significant  at  lower  initial  pressures.  The  same 
behaviour  can  be  seen  in  imperfect  nitrogen  at  much 
higher  shock  Mach  numbers.  The  first  undulation  of 
the  density  ratio  across  a  normal  shock  wave  occurs 
at  shock  Mach  numbers  of  about  10.0  [Ref.  33]  in 
imperfect  air.  This  value  depends  on  the  initial 
pressure.  The  same  undulation,  however,  takes  place 
at  shock  Mach  numbers  of  about  19.0  in  imperfect 
nitrogen  [Ref.  34] . 

It  must  be  noted  that  the  steady  flow  incident 
Mach  number  is  given  by  the  relation  Mo  =  Ms/sin4>0  = 
Ms/cos0^  [Fig.  1].  Consequently,  Mo  becomes  stronger 
with  increasing  0^  for  fixed  Ms.  However,  the  normal 
flow  Mach  number  to  the  incident  shock  wave  is  the 
same  as  M$.  Therefore,  the  flow  state  in  the  region 
behind  the  incident  shock  wave  (region  1)  is  the 
same  for  a  given  Ms  for  every  type  of  reflection  as 
shown  in  Table  1 .  It  is  clear  from  Table  1  that 
the  flow  states  in  the  regions  behind  the  reflected 
shock  wave  (region  2)  and  the  Mach  stem  (region  3) 
depend  on  the  type  of  reflection.  It  was  estimated 
that  oxygen  and  nitrogen  begin  to  dissociate  at 
Ms  =  6.1  and  Ms  =  12.8,  respectively,  in  region  1 
at  P0  =  15  torr  and  T0  =  300  K.  In  region  2  behind 
the  reflected  shock  wave,  dissociations  of  oxygen 
and  nitrogen  occur  at  Ms  =  5.0  and  Ms  =  10.0,  res¬ 
pectively,  along  line  B  (RR  £  MR  transition  line  for 
imperfect  air).  On  the  other  hand,  the  dissocia¬ 
tions  begin  at  Ms  =  6.0  and  Ms  =  12.6,  respectively, 
along  line  D  (CMR  t  DMR  transition  line  for  imperfect 
air).  Ionization  of  oxygen  atoms  and  nitrogen  atoms 
(we  assumed  the  same  ionization  energy  for  0  and  N) 
does  not  occur  until  Ms  =  20.0  in  region  1,  but  it 
takes  place  at  Ms  =  19  5  in  region  2  only  along  line 
B. 

It  should  be  emphasized  that  for  a  reflected 
shock  wave  R  to  occur,  the  flow  behind  the  incident 
shock  wave  I  must  be  supersonic  with  respect  to  the 
reflection  point  P  (in  the  case  of  RR) ,  or  the  first 
triple-point  T  (in  the  case  of  MR) .  Lines  G  and  H 
are  the  lines  corresponding  to  Mi  =1.0  for  perfect 
and  imperfect  air,  respectively.  Under  these  lines, 
i.e.,  in  the  NR  domain,  shock-wave  reflection  does 
not  exist.  It  is  reasonable  to  consider  that  these 
lines  coincide  with  the  lines  9w  =  0,  where  the  first- 
triple-point  trajectory  angle  x  has  its  maximum 
value  Xmax  for  a  fixed  Ms .  If  Ms  is  fixed,  x  is  a 
decreasing  function  with  increasing  8^  (or  9W),  and 
approaches  zero  at  9^  =  9W  =  0wrr-mR  (the  critical 
wedge  angle  for  the  RR  £  MR  transition).  Conse¬ 
quently,  the  NR  domain  does  not  exist  in  practice. 

That  is,  MR  occurs  even  when  the  value  of  8W  is 
very  close  to  zero.  We  will  explain  this  later  in 
detail  in  relation  to  the  variation  of  x  with  9w- 
Lines  E  and  F  represent  the  SMR  £  CMR  transition 
boundary  lines  for  perfect  and  imperfect  air,  res¬ 
pectively.  The  imperfect  line  F  approaches  the 
imperfect  NR  line  H  (where  Mi  =  1.0)  for  Ms  >  10.0. 
This  means  that  the  flow  Mach  number  Mi  behind  the 
incident  shock  wave  decreases  to  nearly  1.0  (but 
must  be  greater  than  1.0),  and  the  reflected  shock 
wave  becomes  a  very  weak  shock  wave.  For  instance, 

Mi  =  1.039  at  Ms  =  10.0  and  Mi  =  1.023  at  Ms  =  16.0 


along  the  line  F. 

It  is  convenient  for  experimental  purposes  to 
draw  the  boundaries  and  domains  corresponding  to 
the  various  types  of  reflection  in  terms  of  the 
incident  shock  Mach  number  Ms,  and  the  actual 
wedge  angle  0W,  i.e.,  in  the  (Ms-6w)-plane.  The 
formation  and  termination  criteria  for  RR,  SMR, 

CMR  and  DMR  in  the  (MS-0W) -plane  are  shown  in  Figs. 
4(a)  and  4(b).  The  solid  lines  are  for  perfect 
air  with  y  =  1.40,  while  the  dashed  lines  are  for 
imperfect  air  at  P0  =  15  torr  and  To  =  300  K.  As 
might  be  expected,  every  boundary  line  for  imperfeci 
air  is  shifted  down  from  its  perfect-gas  value  in 
the  (Ms-0w)-plane  owing  to  real-gas  effects.  This 
result  is  analogous  to  the  results  for  gases  with 
lower  values  of  y  [cf.  Figs.  5,  6].  The  same  type 
of  undulations  of  the  imperfect  gas  lines  occur 
in  the  (MS-0W) -plane  as  in  the  (Ms-0j,) -plane .  The 
divergence  between  the  perfect  and  imperfect  boun¬ 
dary  lines  becomes  greater  with  increasing  M;.  For 
example,  at  Ms  =  5.0,  the  corresponding  actual 
wedge  angle  9W  to  the  CMR  £  DMR  transition  boundary 
is  21.1°  for  perfect  air  and  is  24.7°  for  imperfect 
air  or  22.1°  for  perfect  air  and  13.0°  for  imper¬ 
fect  air  at  Ms  =  10.0.  At  higher  Mach  numbers 
(especially  beyond  Ms  =  10.0),  the  perfect-gas 
lines  level  out  so  that  they  become  nearly  inde¬ 
pendent  of  Ms.  This  can  best  be  understood  by 
referring  to  the  shock  polar  in  the  (4»-9) -plane 
[see  Liepmann  and  Roshko,  Ref.  35],  It  is  seen 
that  for  Mo  "  at  detachment  the  deflection 
angle  has  a  value  9m  =  45.6°  and  the  wave  angle 
^6m  =  67.7°  for  y  =  1.4.  For  Mo  >  10  these  values 
remain  almost  the  same.  However,  for  a  reflected 
wave,  M2  >  1,  then  a  new  polar  has  to  be  drawn  for 
state  2,  which  must  have  the  same  deflection  or 
01  =  02-  Consequently,  0w  =  90°  -  $0  will  achieve 
a  limiting  value  when  0W  't  50° at  Ms  >  10.  Under 
these  conditions,  Ms  =  6.43,  Mq  =  10.0,  0j  =  02  = 
30.7°,  $0  =  40°,  M2  =  2.6  and  <J>j  =  65°.  The  imper¬ 
fect-gas  lines  are  strongly  dependent  on  Ms  due  to 
real -gas  effects.  Every  imperfect-gas  boundary 
line  in  the  (MS-9W) -plane,  as  shown  in  Figs.  4(a) 
and  4(b),  begins  to  undulate  at  high  Mach  number 
but  they  are  always  located  below  the  corresponding 
perfect -gas  boundary  line. 

It  is  worth  noting  that  the  SMR  £  CMR  transi¬ 
tion  line  (M2t  =  1.0)  for  imperfect  nitrogen  at 
Po  =  15  torr  and  T0  =  300  K  will  be  found  in  the 
(MS-0W) -plane  between  the  corresponding  line  for 
a  perfect  diatomic  gas  with  y  =  1.40  and  the  line 
for  imperfect  air  at  Po  =  15  torr  and  To  =  300  K. 
However,  the  SMR  £  CMR  and  CMR  £  DMR  transition 
lines  for  imperfect  nitrogen  at  Po  =  15  torr  and 
To  =  300  K  obtained  by  Ben-Dor  [Refs.  1,  2]  were 
drawn  incorrectly  [see  also  Ref.  36] . 

In  order  to  clarify  the  effects  of  real  gases 
on  the  oblique  shock -wave-reflection  configurations 
it  is  instructive  to  examine  the  following  perfect- 
gas  cases:  (1)  monatomic  gas  (y  =  1.667),  (2)  di¬ 
atomic  gas  (y  =  1.40),  (3)  triatomic  gas  (C02  with 
y  =  1.290),  and  (4)  polyatomic  gas  (SF(,  with  y  = 
1.093).  It  is  well  known  [Ref.  37]  that  an  analogy 
can  be  made  between  real  gases  and  perfect  gases 
with  lower  values  of  y.  It  will  be  seen  that  the 
boundary  lines  for  real  gases  and  those  for  lower 
y  exhibit  significant  similarities. 

Figure  5  shows  the  RR  £  MR  transition  lines 
and  the  SMR  £  CMR  transition  lines  in  the  (Ms-9w)- 
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plane  for  perfect  gases  with  different  y ,  and 
imperfect  air.  The  CMR  £  DMR  transition  lines  are 
shown  in  Fig.  6.  It  is  evident  from  these  figures 
that  each  transition  line  for  a  perfect  gas  drops 
significantly  with  decreasing  y.  This  can  also  be 
understood  by  referring  to  the  shock  polar  in  the 
(<t>-0) -plane  [Ref.  35].  If  y  is  reduced  then  the 
deflection  angle  0m  increases  and  the  wave  angle 
‘t’Sjn  also  rises  until,  at  y  =  1 ,  hypothetically  the 
shock  wave  lies  on  the  body  and  @m  =  $0m  =  00°. 
However,  this  time  the  limiting  M2  ■*  0.  Conse¬ 
quently,  as  y  is  reduced  in  value  0K  is  also  re¬ 
duced  in  order  to  meet  the  condition  6i  =  02  for 
the  lowered  value  of  M2.  This  effect  is  clearly 
illustrated  in  Figs.  5  and  6.  All  perfect-gas  lines 
level  out  at  higher  Mach  numbers  resulting  in  a 
situation  in  which  the  perfect-gas  lines  are  nearly 
independent  of  Ms .  As  Ms  increases,  the  specific- 
heat  ratio  in  imperfect  air  decreases  due  to  the 
vibrational  excitation  and  the  imperfect  lines, 
as  a  result,  begin  to  drop  from  their  perfect-gas 
values  with  y  =  1.40.  At  higher  Ms  dissociation 
of  oxygen  and  nitrogen  proceeds  and  the  specific- 
heat  ratio  (and  the  equilibrium  ye)  of  imperfect 
air  decreases  further.  The  behaviour  of  the  imper¬ 
fect  air  lines  reveals  this  tendency  very  well  as 
can  be  seen  in  Figs.  5  and  6.  Consequently,  the 
lines  are  shifted  downward  with  increasing  Ms  and 
decreasing  y. 

It  is  worthwhile  to  examine  in  detail  the  first 
triple-point  trajectory  angle  y  (in  the  cases  of 
SMR,  CMR  and  DMR),  and  the  second  triple-point- 
trajectory  angle  y*  (in  the  case  of  DMR) .  We  have 
used  the  best  available  analytical  methods  to  date 
for  the  prediction  of  y  developed  by  Law  and  Glass 
[Ref.  19]  and  for  obtaining  y'  developed  by  Ben-Dor 
[Ref.  38].  In  their  analyses,  they  assumed  that 
the  Mach  stem  M  was  straight  and  perpendicular  to 
the  wedge  surface,  and  that  the  second  triple-point 
T',  with  respect  to  the  first  triple-point  T,  moved 
with  the  same  horizontal  velocity  as  the  induced 
flow  behind  the  incident  shock  wave  I .  These 
assumptions  were  deduced  empirically  based  on  the 
geometrical  relation  for  the  Mach  stem  and  on  an 
approximation  for  the  relative  motion  of  the  second 
triple-point.  It  must  be  pointed  out,  however, 
that  these  two  assumptions  do  not  hold  over  the 
entire  ranges  of  0W  and  Ms  as  mentioned  in  Chapter 
2. 

The  values  of  y  and  y '  as  a  function  of  Ms,  with 
the  actual  wedge  angle  9W  as  a  parameter,  are  shown 
in  Figs.  7(a)  and  7(b)  up  to  Ms  =  20.0  for  perfect 
air  with  y  =  1.40.  In  the  following  figures,  the 
line  P*  represents  the  line  M2K  =  1-0,  i.e.,  DMR 
does  not  exist  in  the  region  above  the  line  P*. 

The  values  of  y'  in  this  region  were  merely  obtained 
analytically  and  have  no  physical  meaning.  Figures 
7(a)  and  7(b)  show  that  y'  is,  in  general,  a  de¬ 
creasing  function  with  increasing  9W  and  Ms.  How¬ 
ever,  y'  becomes  nearly  independent  of  Ms  at  higher 
Mach  numbers .  The  same  tendency  can  be  seen  for  y 
at  actual  wedge  angles  0w  less  than  30°.  At  wedge 
angles  greater  than  30°,  the  value  of  y  becomes 
almost  constant  over  the  entire  region  of  Ms.  Ap¬ 
parently,  y  and  y'  are  simple  functions  of  the 
actual  wedge  angle  0W  at  higher  Mach  numbers,  i.e., 

X  *  X(6w)  and  X'  =  X'(9w)-  It  can  be  clearly  seen 
that  for  a  given  0W  the  value  of  y'  is  always 
greater  than  that  of  y  for  perfect  air  with  y  =  1.40. 

The  dependence  of  y  and  y'  on  the  incident 


shock  Mach  number  Ms,  for  constant  0W  for  imperfect 
air  at  P0  =  15  torr  and  To  =  300  K,  is  shown  in 
Figs.  8(a)  and  8(b).  At  low  Mach  numbers  the  value 
of  y  and  y'  approach  those  of  perfect  air  with 
y  =  1.40.  The  imperfect -air  values  of  y  and  y' 
decrease  significantly  with  increasing  Ms  due  to 
real-gas  effects.  Undulations  are  seen  again  in 
the  values  of  y  and  y'  at  higher  shock  Mach  numbers. 
It  is  evident  that  the  perfect-air  values  of  y  and 
y'  are  greater  than  the  imperfect-air  values.  The 
value  of  y'  decreases  much  more  rapidly  than  that 
of  y  as  Ms  increases,  and  becomes  smaller  than 
the  value  of  y  at  higher  Ms  in  the  case  of  imper¬ 
fect  air.  As  noted  above  the  line  P*  represents 
where  M2y  =  1  and  along  the  line  Q*,  uj'  =  0,  as 
explained  below. 

The  reflection  angle  o'  is  defined  as  an  angle 
between  the  wall  surface  and  the  reflected  shock 
wave  in  the  case  of  RR  [Fig.  1(a)],  and  as  an  angle 
between  the  trajectory  path  of  the  first  triple¬ 
point  and  the  reflected  shock  wave  in  the  case  of 
MR  [Figs.  1(a),  2(a),  2(b)].  The  reflection  angle 
o'  is  given  as: 

w  =  (9) 

where  $1  is  the  incident  flow  angle  to  the  reflected 
shock  wave,  and  01  is  the  flow-deflection  angle 
behind  the  incident  shock  wave.  The  reflection 
angle  o'  is  always  positive  for  RR,  SMR  and  CMR. 

The  shift  of  o'  from  a  positive  to  a  negative  value 
with  increasing  Ms  occurs  only  in  the  DMR  domain 
for  imperfect  air  or  gases  with  low  values  of  y. 
Figures  9(a)  and  9(b)  show  u'  as  a  function  of  Ms 
with  0w  as  a  parameter  for  perfect  air  with  y  =  1.40. 
It  is  evident  from  these  figures  that  w'  is  a  de¬ 
creasing  function  of  Ms  at  lower  Mach  numbers  and 
becomes  independent  of  Ms  at  higher  Mach  numbers. 

It  should  be  noted  that  oj'  always  has  a  positive 
value  in  the  case  of  perfect  air  with  y  =  1.40. 
However,  quite  different  features  of  a)'  can  be  seen 
for  imperfect  air.  Figures  10(a)  and  10(b)  show 
the  variation  of  to'  with  Ms  up  to  Ms  =  20.0  for 
imperfect  air.  At  lower  shock  Mach  numbers,  the 
variation  of  w’  is  similar  to  the  perfect-air  case. 
However,  as  Ms  increases,  00 '  decreases  rapidly  and 
becomes  negative  especially  for  larger  0W.  The 
constant  0W  lines  undulate  and  cross  over  owing  to 
real -gas  effects.  These  lead  to  changes  in  the 
incident  and  flow-deflection  angles,  which  are  less 
sensitive  to  real-gas  effects. 

Calculations  were  performed  for  y,  y'  and  w'  as 
a  function  of  Ms  with  0W  as  a  parameter  for  perfect 
gases  with  different  y.  Figure  11  shows  the  varia¬ 
tion  of  y  and  y'  with  Ms  for  a  monatomic  gas  with 
y  =  1.667,  The  same  tendencies  which  appeared  for 
a  diatomic  gas  with  y  =  1.40  can  be  seen  in  Fig.  11. 
The  line  P*  represents  the  line  M2K  =  1.0.  It  can 
be  seen  that  the  value  of  y1  is  always  greater  than 
that  of  y,  i.e.,  the  reflection  angle  w'  is  always 
positive  over  the  entire  range  of  Ms,  as  shown  in 
Fig.  12.  However,  the  values  of  x'  become  less  than 
those  of  y,  and  u>'  has  a  negative  value  above  some 
Mach  numbers  in  the  case  of  a  triatomic  gas  (CO2 
with  y  =  1.290),  as  shown  in  Figs.  13  and  14.  This 
tendency  becomes  more  significant  with  decreasing 
y.  As  can  be  seen  in  Figs.  15  and  16,  x'  becomes 
less  than  x  and  w*  has  a  negative  value  for  a  poly¬ 
atomic  gas  (SFfc  with  y  =  1.093)  at  smaller  Ms  for 
the  same  actual  wedge  angle  6W,  compared  with  the 
triatomic  gas  case.  The  reflection  angle  u'  is 
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shown  in  Fig.  17  as  a  function  of  0W  for  a  given 
Mach  number  Ms  =  10.0.  From  this  figure  it  is 
clear  that  the  reflection  angle  u>'  is  a  very  sensi¬ 
tive  function  of  y  and  a  decrease  in  y  lowers  the 
value  of  j '  significantly  and  shifts  the  value  of 
ui'  into  the  negative  region.  An  increase  in  Ms 
for  a  constant  y  has  the  same  tendency  except  that 
us'  reaches  a  minimum  value  for  a  given  Ms  and  y. 
However,  for  a  perfect  gas  with  y  larger  than  or 
equal  to  1.40,  u>'  cannot  become  negative.  The 
variation  of  ui'  was  made  clear  by  Gvozdeva  et  al 
[Ref.  16] .  They  calculated  the  reflection  angle 
a)'  as  a  function  of  the  incident  shock  velocity 
U0,  and  the  angle  of  incidence  4>0  for  several 
perfect  gases  with  different  y,  and  for  vibration- 
ally  excited  C02  and  air. 

Here  we  may  classify  DMR  into  two  types  and 
define  these  as  follows:  DMR  with  positive  ai*,  as 
illustrated  in  Fig.  18(a),  is  defined  as  a  DMR+  and 
DMR  with  negative  ai',  as  illustrated  in  Fig.  18(c), 
is  defined  as  DMR-.  In  DMR"  the  second  triple-point 
T'  is  located  below  the  first  triple-point  trajec¬ 
tory  path.  The  case  where  X  =  X'  or  u'  =0  appears 
in  Fig.  18(b).  Figure  18(d)  would  represent  a  hypo¬ 
thetical  case  where  oj'  <  0  and  x'  =  0-  The  cal¬ 
culated  results  in  the  present  study  may  be  summar¬ 
ized  as  follows:  DMR+  occurs  all  the  time  for  per¬ 
fect  gases  for  y  £  1.40,  i.e.,  in  a  perfect  diatomic 
gas  and  a  perfect  monatomic  gas.  On  the  other  hand, 
the  transition  from  DMR+  to  DMR"  takes  place  with 
increasing  Ms  in  a  perfect  triatomic  gas  and  in  a 
perfect  polyatomic  gas.  The  same  transition  can  be 
observed  in  the  case  of  an  imperfect  diatomic  gas 
owing  to  real -gas  effects  at  high  Ms. 

Table  2  shows  the  shock  Mach  number  Ms  corre¬ 
sponding  to  in’  =  0  [Fig.  18(b)]  with  0W  as  a  para¬ 
meter.  It  is  clear  from  Table  2  that  the  smaller 
y  gives  a  smaller  Ms  corresponding  to  w'  =0  for  a 
given  actual  wedge  angle  0W.  Line  Q*  in  Figs.  8(a), 
13  and  15  represents  the  line  of  w'  =  0  in  the 
(Ms  "  X.  X')~planes.  The  lines  of  w'  =  0,  i.e., 
the  transition  lines  between  DMR+  and  DMR”  are  shown 
in  the  (Ms-0w) -plane  in  Fig.  19  for  perfect  CO2, 
perfect  SF6  and  imperfect  air  at  Po  =  15  torr  and 
To  =  300  K.  It  must  be  emphasized  that  this 
DMR+  2  DMR"  transition  line  does  not  appear  in  the 
(Ms~0w) -plane  or  the  (Ms-@w) -plane  for  perfect  air. 
The  line  R*  in  Figs.  3(a)  and  4(a)  represents  the 
DMR+  *  DMR"  transition  line  for  imperfect  air.  This 
is  an  important  fact  that  may  be  valuable  in  esti¬ 
mating  the  contribution  of  real -gas  effects  to  the 
shock-wave  configuration  in  the  case  of  air.  The 
DMR+  *  DMR"  transition  occurs  at  decreasingly 
lower  shock  Mach  numbers  as  y  decreases  [see 
Fig.  19].  Ben-Dor  [Ref.  38]  developed  an  analytic¬ 
al  method  for  the  prediction  of  the  second- 
triple-point-trajectory  angle  x'»  and  calculated  x' 
for  perfect  nitrogen  with  y  =  1 .40  and  imperfect 
nitrogen.  Unfortunately,  he  made  a  mistake  in 
explaining  his  result  for  x'-  The  solid  lines  in 
Figs.  7,  8,  9  and  11  of  Ref.  38  must  represent  the 
results  for  perfect  nitrogen,  contrary  to  his  ex¬ 
planation,  and  the  dashed  lines  correspond  to  im¬ 
perfect  nitrogen  at  P0  =  15  torr  and  To  =  300  K. 

The  discussion  which  then  ensues  was  based  on  this 
error,  and  the  comparison  between  the  values  of  x' 
and  x  was  not  made.  However,  we  can  estimate  Ms 
corresponding  to  x1  =  X  ( i - e - ,  w'  =  0),  for  a  given 
0W,  in  the  case  of  imperfect  nitrogen  by  comparing 
the  values  of  x  in  Fig'  31  of  Ref.  1  with  the  values 
of  x'  in  Fig.  9  of  Ref.  38.  The  line  of  w'  =  0  in 
the  (Ms-0w) -plane  for  imperfect  nitrogen  at  PQ  =  15 


torr  and  T0  =  300  K,  as  shown  in  Fig.  19,  is  located 
in  a  much  higher  Ms  region  for  constant  •>« ,  compared 
with  the  present  imperfect-air  result. 

From  the  facts  described  above,  a  diatomic  gas, 
especially  air,  may  be  considered  to  be  an  appro¬ 
priate  test  gas  for  estimating  real-gas  effects  on 
the  shock-wave  configurations.  If  we  can  observe  a 
DMR"  experimentally  at  high  Ms  (for  Ms  >  6  at  P0  = 

15  torr  and  T0  =  300  K) ,  which  does  not  appear  at 
all  in  the  case  of  a  perfect  diatomic  gas,  it  con¬ 
firms  real-gas  effects  on  the  shock-wave  configur¬ 
ations  . 

As  may  be  seen  in  Fig.  15,  the  second  triple¬ 
point-trajectory  angle  x'  becomes  zero  at  Gw  lar¬ 
ger  than  40°  for  Ms  <  10  in  the  case  of  perfect 
SFfe  with  y  =  1.093.  This  means,  as  illustrated  in 
Fig.  18(d),  that  the  second  triple-point  T'  would 
attach  to  the  wedge  surface,  i.e.,  the  reflected 
shock  wave  R  would  strike  the  wedge  surface  and  be 
reflected  like  a  RR.  The  second  triple-point  T’ 
attaches  to  the  wedge  surface  at  Ms  =  6.5  for 
6W  =  40°  and  at  Ms  =  5.0  for  0W  =  45°.  However,  it 
should  be  pointed  out  that  these  points  are  above 
the  MR  termination  line  in  the  (MS-9W) -plane  for 
perfect  SF^  as  can  be  seen  in  Fig.  5.  (®wRR-MR  *s 
40.1°  at  Ms  =  5.0  and  9wrr-mr  is  38.9°  at  Ms  =  6.5, 
where  9wrr_mr  is  a  critical  angle  for  the  RR  £  MR 
transition.)  Therefore,  a  reflection  of  the  type 
shown  in  Fig.  18(d)  cannot  exist  in  practice.  This 
will  be  explained  subsequently  in  relation  to  the 
two-  and  three-shock  theories. 

As  noted  previously,  there  are  mainly  two 
different  criteria  which  exist  for  the  termination 
of  RR.  The  first  one  is  known  as  the  detachment 
criterion,  and  the  second  one  is  referred  to  as 
the  mechanical  equilibrium  criterion  which  was 
found  experimentally  to  be  inapplicable  to  nonsta¬ 
tionary  flows.  Note  that  the  detachment  criterion 
is  based  on  the  two-shock  theory  [Ref.  5]  and  the 
mechanical  equilibrium  criterion  originates  in  the 
three-shock  theory  [Ref.  6].  In  order  to  clarify 
the  difference  between  these  two  criteria,  the 
termination  lines  of  RR(MR)  in  the  (Ms-0w)-plane 
are  drawn  as  a  function  of  Ms  in  Fig.  20  for  a 
perfect  diatomic  gas  with  y  =  1.40.  As  is  evident 
from  Fig.  20,  the  actual  wedge  angle  9W  correspond¬ 
ing  to  the  mechanical  equilibrium  criterion  is 
larger  than  the  detachment  criterion  over  the  entire 
range  of  Ms .  In  the  region  between  the  detachment 
and  mechanical  equilibrium  criteria,  MR  may  exist 
if  we  apply  the  mechanical  equilibrium  criterion 
as  the  MR  termination  criterion.  In  other  words, 
we  can  obtain  analytically  MR  with  a  finite  value 
of  x  in  this  region  using  the  three-shock  theory. 
However,  the  detachment  criterion  must  be  applied 
as  the  MR  termination  criterion  for  nonstationary 
flows.  Therefore,  MR  cannot  exist  above  the  MR 
termination  line  determined  by  the  detachment 
criterion  for  the  nonstationary  case.  Consequently, 
if  Ms  is  fixed  while  increasing  the  actual  wedge 
angle  0W,  the  first  triple-point-trajectory  angle  x 
decreases  monotonically  and  must  drop  to  zero 
discontinuously  at  0W  =  OwRR-MR-  Note  that  x  and 
X'  were  calculated  based  on  three-shock  theory. 

The  points  corresponding  to  x'  =0  are  above  the  MR 
termination  line  in  the  (Ms-9w)-plane  [Fig.  5]  and 
therefore  they  do  not  exist  in  practice. 

On  the  other  hand,  if  one  starts  at  a  fixed 
shock  Mach  number  Ms  and  decreases  the  actual  wedge 
angle  9w,  then  the  first  triple-point-trajectory 
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angle  x  increases  monotonically,  as  can  be  seen  in 
Figs.  7(a),  8(a),  11,  13  and  15.  The  values  of  x 
as  a  function  of  Ms  for  very  small  0W  are  shown  in 
Fig.  21*  for  a  perfect  diatomic  gas  with  y  =  1.40. 

For  fixed  Ms,  X  increases  monotonically  as  9w 
decreases.  A  dot-dashed  line  in  Fig.  21  represents 
the  line  of  Mj  =  1.0,  which  corresponds  to  the 
critical  line  for  the  existence  of  a  reflected 
shock  wave  R  and  also  may  be  considered  as  the  line 
for  which  the  critical  wedge  angle  9*  =  0°.  Conse¬ 
quently,  the  line  of  x  for  a  given  6W  terminates 
at  the  Mj  =  1.0  line  (the  critical  value  of  x  for 
9W  =  0°),and  when  9W  -*■  0°  this  line  must  merge  with 
the  Mi  =  1.0  line.  As  is  evident  in  Fig.  21,  the 
value  of  x  predicted  by  the  method  of  Law  and  Glass 
(Ref.  19]  monotonically  approaches  the  Mi  =  1.0  line 
with  decreasing  6 w  for  Ms  >  2.0.  However,  the  value 
of  x  becomes  larger  than  the  Ml  =  1.0  line  at  Ms 
smaller  than  2.0  for  very  small  actual  wedge  angles 
(0W  £  5°).  Ben-Dor  [Ref.  1]  proposed  an  alternative 
method  to  predict  x  for  very  small  wedge  angle 
(9W  5  4°)  empirically  as  follows: 

x  =  Vi.o '  ew  (10> 

where  ®wMi=l.o  is  the  effective  wedge  angle  corre¬ 
sponding  to  Mi  =  1.0.  This  method,  however,  is  quite 
an  artificial  one  and  is  inconsistent  with  the  method 
developed  by  Law  and  Glass  [Ref.  19]  for  predicting  x- 
A  better  method  for  predicting  the  correct  value  of 
X  at  smaller  Ms  and  0W  is  still  required. 

3.2  Comparison  of  Analytical  Results  with  Available 
Experimental  Results 

In  the  following,  the  present  analytical  re¬ 
sults  of  pseudo-stationary  oblique  shock-wave 
reflection  are  compared  with  experiment.  Figure  22 
is  a  reproduction  of  Fig.  3(a)  and  represents  the 
domains  and  transition  boundaries  for  various  types 
of  reflection  in  the  (Ms-0,i) -plane.  The  solid 
lines  are  perfect  air  with  y  =  7/5,  while  the 
dashed  lines  are  for  imperfect  air  with  a  given 
initial  pressure  Po  =  15  torr  and  a  temperature 
To  *  300  K.  The  experimental  data  of  Deschambault 
for  air  [Ref.  39]  and  the  data  of  Ben-Dor  for  nitro¬ 
gen  [Ref.  40]  are  plotted  in  Fig.  22  to  check  the 
present  analysis  for  the  transition  boundaries. 
Deschambault  [Ref.  39]  conducted  interferometric 
experiments  for  air  in  the  UTIAS  hypervelocity 
shock  tube  with  a  10cm  x  18cm  cross-section.  The 
initial  temperatures  T0  were  nearly  300  K  and 
pressures  were  in  the  range  5.0  torr  £  Po  5  785.5 
torr.  The  incident  shock  Mach  number  range  covered 
was  1.2  £  Ms  £  10.3  over  a  series  of  wedge  angles 
2.0°  £  0W  £  63.4°,  The  initial  conditions  and  the 
experimental  results  are  summarized  in  Table  3.  It 
is  evident  from  Fig.  22  that  the  experimental  data, 
except  for  RR  and  one  or  two  other  points,  lie 
inside  the  domains  predicted  by  the  perfect-gas 
model.  It  can  be  seen  that  RR  persists  below  the 
perfect  RR  t  MR  transition  line  determined  by  the 
detachment  criterion.  This  persistence  of  RR  at 
low  shock  Mach  numbers  was  called  "von  Neumann's 
paradox"  [Refs.  4,  9,  11].  It  is  clear  from  the 
experimental  results  for  air  that "von  Neumann's 
paradox"  also  exists  at  higher  shock  Mach  numbers. 

The  possibility  exists  that  RR  occurs  below  the 
analytical  terminating  boundary  line  for  perfect 


air  for  all  Ms.  However,  a  detailed  and  careful 
set  of  experiments  must  be  performed  at  higher  Ms 
for  8W  <  50°  to  verify  this  point.  A  comparison 
of  the  experimental  results  show  good  agreement 
with  all  the  other  perfect-gas  boundary  lines, 
except  for  the  perfect  RR  t  MR  transition  line. 

It  should  be  stressed  that  these  results  for  6^,  do 
not  require  the  measurement  of  x  and  can  be  con¬ 
sidered  more  accurate  than  plots  of  6W  [cf.  Fig. 

23], 

The  experimental  results  of  Deschambault  [Ref. 
39]  for  air  (see  Table  3  for  details)  as  well  as 
the  data  of  Ben-Dor  [Ref.  40]  for  nitrogen  are 
plotted  in  the  (Ms-9w) -plane  in  Fig.  23  [a  repro¬ 
duction  of  Fig.  4(a)],  It  is  noted  from  Fig.  23 
that  several  SMR  points  lie  in  the  domain  corre¬ 
sponding  to  CMR.  For  air,  SMR  were  observed  at 
Ms  =  3.0  for  0W  =  20°  and  at  Ms  =  6.41  for  6w  =  10°. 
Meanwhile,  Ben-Dor  reported  SMR  for  0w  =  10°  at 
Ms  =  5.92,  6.79  and  7.58.  On  the  other  hand, 
several  of  the  data  points  for  air  lie  below  the 
perfect  CMR  *■  DMR  transition  line.  That  is, 
Deschambault  observed  DMR  configurations  at 
Ms  =  3.79  for  0w  =  30°,  at  Ms  =  3.84  and  4.77  for 
8W  =  27°,  and  at  Ms  =  8.85  for  0w  =  20°.  It  was 
not  possible  to  redo  some  of  these  experiments  to 
check  their  repeatability  at  high  Ms  owing  to 
possible  damage  to  the  interferometric-quality 
windows .  Some  experimental  points  for  air  even 
lie  below  the  imperfect  CMR  *  DMR  transition  line. 
These  facts  suggest  that  the  existing  criteria  for 
the  SMR  ^  CMR  transition  (M2t  =  1.0)  and  the 
CMR  *  DMR  transition  (M2g  =  1.0)  might  be  neces¬ 
sary  but  not  sufficient  conditions  for  these 
transitions  [see  Ref.  36] .  Further  investigation 
regarding  this  point  is  necessary.  Apart  from 
this  disagreement  and  the  persistence  of  RR  below 
the  detachment  criterion,  the  experimental  points, 
in  general,  lie  inside  of  their  domains  predicted 
by  the  perfect-gas  model  in  the  (MS-0W) -plane  [Fig. 
23].  Consequently,  the  analysis  as  given  in  Refs. 
1-3  is  reasonable,  if  not  precise,  in  modelling 
the  actual  physics  of  the  problem. 

Figure  24  shows  the  experimental  data  in  the 
(Ms-9w) -plane  for  air,  nitrogen  and  oxygen  from 
several  sources  [Refs.  9,  18,  19,  39  and  40]  in 
the  range  1.0  <  Ms  £  3.0.  It  is  evident  from  the 
data  of  Smith  [Ref.  9]  that  RR  persists  even  below 
the  imperfect  detachment  line  at  P0  =  15  torr  and 
T0  =  300  K.  Note  that  the  contribution  of  real- 
gas  effects  is  merely  due  to  vibrational  excita¬ 
tion  of  oxygen  and  nitrogen  molecules  in  this  shock 
Mach  number  range.  Although  Smith  had  noticed  a 
CMR  configuration,  he  did  not  distinguish  between 
SMR  and  CMR,  and  referred  to  them  as  MR.  Except 
for  RR,  the  experimental  points  appear  to  lie  well 
inside  of  the  analytically  predicted  domains  for 
Ms  £  3.0.  However,  had  Smith  distinguished  between 
SMR,  CMR  and  DMR  we  might  have  seen  some  discrep¬ 
ancies.  For  example,  Deschambault  [Ref.  39] 
observed  a  DMR  configuration  at  Ms  =  2.08  for 
0W  =  49°,  in  which  the  data  point  was  located 
in  the  CMR  domain.  This  fact  supports  the  hypo¬ 
thesis  of  Ando  and  Glass  [Ref.  20]  for  results  in 
C02  that  the  SMR  t  CMR  and  CMR  t  DMR  transition 
lines  meet  at  the  RR  t  MR  transition  boundary. 
However,  more  experimental  data  in  the  CMR  domain 
for  9w  £  40°  are  required  in  order  to  have  addi¬ 
tional  evidence. 

Figure  25  shows  the  predicted  values  of  x  and 
X'  versus  Ms  for  0w  =  27°  and  40°  and  their 
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comparison  with  experiment.  In  this  figure  the 
solid  lines  are  for  perfect  air  with  y  =  7/5, 
while  the  dashed  lines  are  for  imperfect  air  at 
Po  =  15  torr  and  To  =  300  K.  The  experimental 
results  for  air  from  Ref.  39  and  nitrogen  from 
Ref.  40  are  plotted.  It  is  clear  from  Fig.  25 
that  the  experimental  data  of  x  and  x'  for  9W  = 

27°  agree  reasonably  well  with  the  perfect-gas 
values  over  the  entire  range  of  Ms.  However,  for 
Ok  =  40°  the  agreement  is  only  good  for  Ms  <  6 . 

For  the  range  Ms  2  6  the  agreement  becomes  worse, 
and  the  experimental  values  of  x  and  x'  deviate 
from  the  values  predicted  by  the  perfect-gas  model . 

It  is  worth  mentioning  that  the  method  for  predict¬ 
ing  x  proposed  by  Law  and  Glass  [Ref.  19]  is  only 
good  in  the  range  0w  <  40°  [Bazhenova  et  al.  Ref. 

18],  and  that  the  values  of  x  and  x*  are  relatively 
small  at  higher  wedge  angles  which  might  lead  to 
noticeable  experimental  error.  The  calculated  and 
measured  values  of  reflection  angle  os'  are  plotted 
as  a  function  of  Ms  for  0W  =  27°  and  40°  in  Fig. 

26.  The  experimental  points  for  air  lie  close  to 
the  solid  lines  (perfect  gas)  with  y  =  7/5.  It 
should  be  noted  that  DMR",  which  may  occur  in  imper¬ 
fect  air,  and  the  DMR+  t  DMR~  transition  (o'  =0) 
depend  on  the  initial  thermodynamic  conditions.  Such 
transitions  were  not  observed  in  experiments  in  air 
[Ref.  39]  or  nitrogen  [Ref.  40]  in  the  range  Ms  £  10. 
That  is,  x'  is  always  greater  than  x.  and  therefore 
oj'  is  always  positive  [see  Table  3  of  present  report 
and  Table  2(a)  of  Ref.  40],  again  confirming  that 
real-gas  effects  were  not  significant. 

As  mentioned  previously,  the  transition  from 
DMR+  to  DMR-  occurs  with  increasing  Ms  in  a  perfect 
gas  with  y  less  than  1.40.  The  calculated  values  of 
reflection  angle  w'  for  a  perfect  gas  with  y  =  1.290 
were  compared  with  the  available  experimental  data 
for  C02  [Table  2  of  Ref.  13]  in  Fig.  27.  It  is 
evident  from  Fig.  27  that  the  experimental  data  for 
the  reflection  angle  w'  decreases  rapidly  and  be¬ 
comes  negative  for  Ms  ;  4.  Ando  plotted  the  differ¬ 
ence  Aw'  between  the  predicted  and  measured  u>'  for 
CO2  [see  Appendix  F  of  Ref.  13] .  He  showed  that  the 
predicted  value  of  id'  based  on  the  perfect-gas  model 
with  y  =  1.290  deviated  from  the  measured  values  sig¬ 
nificantly  for  Mj.  >  3.5.  The  quasi-equilibrium 
models,  taking  into  account  the  excitation  of  the 
vibrational  modes,  agreed  best  with  experiment  for 
<*>'  [Ref.  13],  since  real-gas  effects  contribute  at 
Ms  >  3.5  in  C02- 

This  is  confirmed  in  Ref.  36  where  the  vibra¬ 
tional  effects  on  u>’  are  also  plotted  in  a  differ¬ 
ent  form.  The  experimental  data  is  in  good  agree¬ 
ment  with  this  plot,  as  shown  in  Fig.  27  for  the 
case  of  oj*  for  a  set  of  wedge  angles  of  20°, 

30°  and  40°  for  CO2  in  vibrational  equilibrium. 

Better  agreement  is  obtained  at  the  higher  angles 
since  the  analytical  curves  lie  closer  together  in 
that  range.  In  addition,  the  fact  that  the  plot 
has  0W  rather  than  0W  as  a  variable  also  adds  to 
deviations  since  x  nwst  be  known  to  obtain  8W  but 
not  0^.  Consequently,  any  errors  in  the  analytical 
results  for  x  will  also  show  up  in  this  plot.  On 
the  whole  though,  the  agreement  is  quite  good. 
Consequently,  even  though  the  types  of  reflections 
agree  best  with  y  =  1.290  in  Ref.  13,  the  measure¬ 
ments  of  w'  are  much  more  basic  as  they  agree  with 
the  relaxation  lengths.  This  raises  a  question  as 
to  the  validity  of  the  perfect-gas  transition  lines 
for  CO2  in  the  (Ms-0w)  and  (Ms-0w)  planes. 

It  is  worthwhile  to  examine  the  relaxation 


lengths  behind  the  shock  waves.  The  relaxation 
length  for  vibrational  or  dissociational  excitation 
which  depends  on  the  gas  used  and  the  initial 
thermodynamic  conditions,  is  an  important  factor 
in  estimating  the  state  of  the  gas  behind  the 
shock-front.  If  the  relaxation  length  behind  a 
shock  wave  is  large  compared  with  the  shock  thick-  | 
ness,  which  is  on  the  order  of  a  few  gas  mean-free- 
paths,  the  state  of  the  gas  at  the  shock  front 
can  be  considered  as  frozen.  On  the  other  hand, 
if  the  relaxation  length  is  negligibly  small 
(order  of  mean-free-paths) ,  the  state  of  the  gas 
behind  the  shock-front  is  in  equilibrium.  Various 
reports  exist  dealing  with  vibrational  and  disso¬ 
ciational  relaxations  in  N2,  O2  and  air  [Refs. 

41-44] .  The  relaxation  lengths  behind  a  normal 
shock  wave  versus  Ms  are  plotted  in  Fig.  28  for 
the  various  excitation  modes  (vibrational  excita¬ 
tion  for  O2  and  N2,  and  dissociational  excitation 
for  O2,  N2  and  overall  air).  These  curves  are 
for  the  initial  conditions  of  PQ  =  10  torr  and 
T0  =  300  K,  and  were  obtained  on  the  basis  of  the 
above  theoretical  and  empirical  analyses  [Refs. 

41,  43,  44].  From  Fig.  28,  we  can  estimate  the 
relaxation  lengths  behind  shock  waves  for  the 
experimental  cases  of  air  [Ref.  39]  and  nitrogen 
[Ref.  40],  Table  4  shows  the  vibrational  and 
dissociational  relaxation  lengths  behind  the 
incident  shock  wave  in  the  case  of  air  [Table  3] 
for  Ms  >  5.  It  is  evident  from  Table  4  that  the 
relaxation  lengths  except  for  02  vibration  are 
large.  The  shortest  relaxation  length  for  02 
vibration  is  0.05  mm  (Exp.  #416),  and  ^02vib  of 
other  cases  are  in  the  range  of  0.15  n,  13  mm. 
Consequently,  we  can  estimate  that  the  state  of 
the  gas  at  the  incident  shock-front  is  nearly 
frozen.  Only  the  02  vibrational  mode  may  be 
partially  excited.  The  mole  fraction  of  oxygen 
molecules  in  air  is  about  20%;  therefore,  the 
partial  excitation  of  02  vibration  may  not  sig¬ 
nificantly  influence  the  resulting  thermodynamic 
state  of  air  at  the  shock-front. 

Table  5  shows  the  vibrational  relaxation  length 
X-N-vib  behind  the  incident  shock  wave  for  the  case 
of^nitrogen  [Table  2(a)  of  Ref.  40]  for  Ms  >  5. 

The  relaxation  lengths  for  N2  vibrational  excita¬ 
tion  are  quite  large  (the  shortest  fy^vib  69  m")  • 
The  relaxation  lengths  behind  the  Macn  stem  (in 
the  case  of  MR)  can  be  considered  to  be  of  the 
same  order  of  magnitude  as  those  behind  the  inci¬ 
dent  shock  wave. 

The  relaxation  lengths  behind  the  reflected 
shock  wave  are  much  shorter  than  those  behind  the 
incident  shock  wave  owing  to  the  two-step  compres¬ 
sions  through  the  incident  and  reflected  shock 
waves.  Even  if  we  assume  the  relaxation  lengths 
behind  the  reflected  shock  wave  to  be  one  order  of 
magnitude  less  than  those  behind  the  incident 
shock  wave,  I^vib  behind  the  reflected  shock  wave 
is  still  large  for  N2  experiments  of  Ben-Dor  [Ref. 
[40].  In  Deschambault's  experiments  in  air  [Ref.  39] 
?02  is  still  small  behind  the  reflected  shock  wave 
when  Ms  >  5  and  p0  >  30  torr  approximately  (Table  4). 
From  these  estimates  we  can  conclude  that  the  state 
of  the  gas  at  each  shock-front  is  frozen  for  N2  ex¬ 
periments  [Ref.  40]  and  is  nearly  frozen,  except  for 
O2  vibrational  excitation  in  some  cases,  for  air  ex¬ 
periments  [Ref.  39],  Consequently,  the  shock  shapes 
around  the  confluence  point  are  determined  at  the 
frozen  shock-fronts,  and  therefore  the  experimental 
data  for  air  and  nitrogen  agree,  in  general,  with 
the  analytical  results  based  on  the  perfect-gas  model. 
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The  discrepancy  between  the  perfect-gas  analysis 
and  experiment,  i.e.,  the  persistence  of  RR  below 
the  detachment  criterion  and  the  disagreement  of 
some  experimental  points  with  the  SMR  t  CMR  and 
CMR  t  DMR  transitions,  are  probably  due  to  the 
criteria  themselves  rather  than  to  real -gas  effects 
suggested  earlier  in  Refs.  1-3,  or  possibly  dis¬ 
placement  effects  of  the  boundary  layer  on  the  wall 
as  suggested  by  Homung  et  al  [Ref.  7].  However, 
the  latter  effect  has  hardly  been  verified  and  is 
considered  in  some  detail  in  Ref.  36. 


4.  CONCLUSIONS 


with  experimental  data  for  persistence  of  RR 
and  the  measured  reflected  wave  angle  u)1. 

4.  To  resolve  some  of  the  remaining  problems  such 
as:  a  better  analytical  method  for  predicting 
the  first  triple-point  trajectory  angle  x  for 
very  small  wedge  angle  (9W  5  5°),  and  the 
location  of  the  kink  K  (or  the  second  triple 
point  T')  for  large  wedge  angle  (9w  £  40°). 

Some  of  these  problems  have  already  been  considered 
in  some  detail  and  will  be  reported  on  in  Ref.  36, 
in  the  very  near  future. 


An  analytical  and  numerical  study  was  made  of 
oblique-shock-wave  reflection  in  air  in  pseudo¬ 
stationary  flow.  The  domains  and  transition 
boundaries  between  the  various  types  of  reflection 
(RR,  SMR,  CMR  and  DMR)  were  established  using  the 
transition  criteria  in  R'efs.  1-3,  up  to  Ms  =  20 
for  both  perfect  and  imperfect  air.  Additional 
calculations  were  made  for  perfect  gases  with 
differing  y  in  order  to  clarify  the  effect  of  the 
specific-heat  ratio  y  on  the  shock-wave  configura¬ 
tions.  The  transition  boundary  lines  significantly 
shift  downwards  with  decreasing  y,  as  for  real 
gases.  It  was  verified  that  the  reflection  angle 
oj'  was  a  very  sensitive  function  of  y.  A  decrease 
in  y  lowered  the  value  of  to '  significantly.  A 
decrease  in  y  shifted  the  value  of  to'  towards 
negative  values  under  certain  conditions.  For 
a  perfect  gas  with  y  =  1.67  or  1.40,  it  was  not 
possible  to  have  a  negative  value  of  u)' .  It  was 
found  that  ui'  could  be  positive  or  negative  for 
imperfect  air  at  higher  Ms. 

Transition  boundaries  predicted  on  the  basis 
of  a  perfect  gas  with  y  =  1.40  agree,  in  general, 
with  the  available  experimental  data  for  air  and 
nitrogen.  However,  RR  persists  below  the  perfect 
termination  line  determined  by  the  detachment 
criterion,  and  SMR  and  DMR  can  sometimes  occur 
outside  their  analytically  predicted  domains.  The 
experimental  data  used  in  the  present  study  was 
checked  in  relation  to  the  relaxation  lengths 
behind  the  shock  waves,  which  are  good  indicators 
for  estimating  the  states  of  the  gas  immediately 
behind  the  shock  waves.  It  was  clarified  that  the 
flow  states  at  the  shock  fronts  were  frozen  or 
nearly  frozen  regarding  vibrational  excitation  and 
dissociation.  Consequently,  it  can  be  concluded 
that,  within  the  limit  of  the  available  experi¬ 
mental  data,  the  various  wave  systems  are  deter¬ 
mined  at  the  frozen  shock  fronts  and,  in  general, 
the  perfect-gas  analysis  predicts  well  the  shock- 
wave  configurations.  The  discrepancies  between  the 
present  analysis  and  experiment  are  probably  due 
to  the  basic  assumptions  and  the  transition  criteria 
themselves. 

Four  outstanding  problems  remain  for  future 
consideration: 

1.  Improved  criteria  for  all  transition  lines. 

This  could  explain  the  "von  Neumann  paradox" 
and  improve  the  agreement  between  analysis 
and  experiment. 

2.  To  determine  experimentally  the  persistence  of 
RR  over  a  significant  range  of  3  £  Ms  S  10. 

3.  A  detailed  analysis  of  the  boundary-layer 
growth  on  the  wedge  surface  and  a  comparison 
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Table  1 


Degree  of  dissociation  in  states  (1)  to  (3)  for  oblique  shock-wave  reflections 


for  imperfect  air  at  PQ  =  15  torr  and  To  =  300  K 
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Table  2 


Shock  Mach  number  corresponding  to  U3'=0  for  C02>  SF6  and  imperfect  air 
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Table  3 


Initial  conditions  and  experimental  results  in  air 
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Table  3  -  Continued 


Initial  conditions  and  experimental  results  in  air 
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Table  4 


Relaxation  lengths  (H*)  in  mm  behind  incident  shock  waves  in  air 
for  given  initial  pressures  pQ  at  300  K 
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FIG.  2  SCHEMATIC  DIAGRAM  OF  COMPLEX  AND  DOUBLE  MACH  REFLECTIONS: 

(a)  COMPLEX  MACH  REFLECTION  ( CMR ) ,  (b)  DOUBLE  MACH  REFLECTION  ( DMR ) . 
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FIG.  4(a)  REGIONS  OF  VARIOUS  TYPES  OF  SHOCK-WAVE  REFLECTION  IN  (MS-8W)  PLANE 
SOLID  LINES  ARE  FOR  PERFECT  AIR  WITH  y  =  1.40;  DASHED  LINES  ARE 
FOR  IMPERFECT  AIR  AT  P0  =  15  TORR  AND  T0  =  300  K. 
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G  6  TRANSITION  LINES  FOR  RR  *  MR  AND  CMR  t  DMR  FOR  PERFECT  GASES 
WITH  DIFFERENT  y  AND  IMPERFECT  AIR. 
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FIG.  7(a)  VARIATION  OF  x  AND  x’  WITH  Ms  FOR  GIVEN  0W  (PERFECT  AIR 
WITH  y  =  1.40). 
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UATION  OF  REFLECTION  ANGLE  to' 
.tPERFECT  AIR  AT  Pn  =  15  TORR  Af 
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Ms 

:IG.  16  VARIATION  OF  REFLECTION  ANGLE  u>'  WITH  M,  FOR  GIVEN  9W 
(POLYATOMIC  GAS  -  SF6  WITH  y  =  1.093). 


V/////////////////////////////////////////////A 


(d)  x'=0,  w'<0 


-  CONTINUED  -  SEVERAL  PATTERNS  OF  DOUBLE  MACH  REFLECTION 
CASE  (d)  IS  HYPOTHETICAL. 


<k 


FIG.  19  TRANSITION  LINES  FOR  DMR+  t  DMR”  (w*  =  0)  FOR  C02, 
SFs,  IMPERFECT  AIR  AND  IMPERFECT  NITROGEN. 
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FIG.  20  DETACHMENT  AND  MECHANICAL  EQUILIBRIUM  CRITERIA  FOR 
MR(RR)  TERMINATION,  v  =  1.40. 


FIG.  22 


COMPARISON  OF  PREDICTED  DOMAINS  OF  VARIOUS  TYPES  OF  REFLECTION 
AND  EXPERIMENTAL  RESULTS  IN  THE  (Ms-9w)  PLANE.  SOLID  LINES  ARE 
FOR  PERFECT  AIR  WITH  y  =  1.40.  DASHED  LINES  ARE  FOR  IMPERFECT 
AIR  AT  P0  =  15  TORR  AND  Tg  =  300  K.  SOLID  SYMBOLS  ARE  AIR  DATA 
FROM  OESCHAMBAULT  [REF.  39].  OPEN  SYMBOLS  ARE  N,  DATA  FROM 
BEN-DOR  [REF.  40].  RR  0  ■,  SMR  A  A,  CMR  0  •,  DMR  9  +. 


FIG.  23  COMPARISON  OF  PREDICTED  DOMAINS  OF  VARIOUS  TYPES  OF  REFLECTION 
AND  THE  EXPERIMENT  IN  THE  (MS-6W)  PLANE.  SOLID  LINES  ARE  FOR 
PERFECT  AIR  WITH  y  =  1.40.  DASHED  LINES  ARE  FOR  IMPERFECT  AIR 
AT  P0  =  15  TORR  AND  T0  =  300  K.  SOLID  SYMBOLS  ARE  AIR  DATA 
FROM  DESCHAMBAULT  [REF.  39].  OPEN  SYMBOLS  ARE  N2  DATA  FROM 
BEN-DOR  [REF.  40].  RR  ,  SMR  AA  ,  CMR  0  §,  DMR  ®  +. 
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FIG.  24  COMPARISON  OF  PREDICTED  DOMAINS  OF  VARIOUS  TYPES  OF  REFLECTION  AND 
THE  EXPERIMENT  IN  THE  (MS-9W)  PLANE  IN  THE  RANGE  1.0  <  Ms  S  3.0. 

SOLID  LINES  -  ARE  FOR  PERFECT  AIR  WITH  y  =  1.40.  DASHED  LINES 

.  ARE  FOR  IMPERFECT  AIR  AT  P0  =  15  TORR  AND  T0  =  300  K.  DASH-DOT 

LINE  -  SHOWS  PERSISTENCE  OF  REGULAR  REFLECTION  RR.  EXPERIMENTAL 

DATA  -  AIR:  SMITH  9  ♦  [REF.  9],  BAZHENOVA  ET  AL  n  [REF.  18], 
DESCHAMBAULT  [REF.  39];  NITROGEN:  BAZHENOVA  ET  AL  A  [REF.  18] 

BEN-DOR  Oav  [REF.  40];  OXYGEN:  LAW  AND  GLASS  A  v  [REF.  19].  RR  • 
SMR  A  A  A  A  ,  CMR  ▼  v  ▼  ,  DMR  ■  B  ,  MR  ♦. 


FIG.  25  EXPERIMENTAL  COMPARISON  OF  x  AND  X'  WITH  Ms  FOR  GIVEN  0W.  SOLID 
.INES  SHOW  PERFECT  AIR  WITH  y  =  1.40.  DASHED  LINES  SHOW  IMPERFECT 
AIR  AT  P0  =  15  TORR  AND  T0  =  300  K.  SOLID  SYMBOLS  ARE  AIR  DATA 
FROM  DESCHAMBAULT  [REF.  39].  OPEN  SYMBOLS  ARE  Nz  DATA  FROM  BEN-DOR 
[REF.  40].  ▼:  X'  (0W  =  27°),  x*  (0W  =  40°),  A:  x  (0W  =  27°), 

•  0:  x  (0W  =  40°). 


FIG.  27  COMPARISON  OF  REFLECTION  ANGLE  w'  WITH  SHOCK  MACH  NUMBER  Ms 

FOR  GIVEN  WEDGE  ANGLE  0W  WITH  EXPERIMENTS  IN  C02  FROM  REF.  13. 

0:  0W  =  20°,  a:  0w  =  30°,  □:  0W  *  40° ,  - y  =  1 .29, . 

VIBRATIONAL  EQUILIBRIUM. 


Relaxation  Path  Lengths, 


N2  Dissociation 
V  (Logan) 


02  Dissociation 
(  (Logan) 


N2  Vibration 
^(Vincenti  a 
\J Kruger) 


02  Vibration\  \  Air  Overall 
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FIG.  28  RELAXATION  LENGTHS  l*  BEHIND  NORMAL  SHOCK  WAVE  VERSUS  Ms 
FOR  INITIAL  CONDITIONS  OF  P0  =  10  TORR,  T0  =  300  K  [AFTER 
REFS.  41,  43,  44].  NOTE  THAT  THE  VARIOUS  DEPEND  ON 
THE  ACTUAL  P  AND  T  BEHIND  SHOCK  WAVE. 
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